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1.    SUMMARY: 

This  Technical  Note  establishes  the  relations  for  determining  the 
flow  rates  of  the  Transonic  Turbine  Test  Rig  at  the  Propulsion  Laboratory 
of  the  Department  of  Aeronautics,  U.  S.  Naval  Postgraduate  School. 

Experimental  results  of  flow  nozzle  calibration  tests  and  of  partial 
range  plenum  labyrinth  seal  leak  tests  are  presented.   In  addition,  formulas 
for  the  estimation  of  shaft  seal  leakage  are  included.   The  equations  that 
must  be  used  to  determine  the  flow  rate  through  the  turbine  are  listed  in 
Table  IV. 

The  contents  of  this  Technical  Note  will  be  incorporated  in  a  forth- 
coming thesis  by  the  author. 


OLS  f( 

2 
A  -  Area  (in  ) 

b  -  heighth  of  labyrinth  </al  intertooth  chamber  (in) 

f 

C  -  constant  in  sharp-/.-/ged  orifice  equation 

D  -  diameter  of  pipe  upstream  of  orifice  (nozzle)  (in.) 

D  -  diameter  of  ori/ice-minimum  diameter  of  nozzle  (in.) 
2  ■*? 

2 
g  -  acceleration  due  to  gravity  (32.]  7  lbm  ft/lbs  sec  ) 

b^  -  differential  pressure  (in  H„0) 

K  -  discharge  coefficient  (dimensionless) 

n  -  number  of  throttlings  in  labyrinth  seal  (dimensionless) 

P,  -  pressure  upstream  of  orifice  or  nozzle  (in.  Hg) 

pQ  -  total  pressure  at  labyrinth  seal  inlet  (psia) 

p-,  -  pressure  upstream  of  orifice  or  nozzle  (psia) 

p  -  static  pressure  at  labyrinth  seal  discharge  (psia) 

R  -  gas  constant  for  air  (53.35  ft-lb/lbm°R) 

R  -  Reynold's  Number  based  on  D2  (dimensionless) 

r  -  overall  pressure  ratio  of  labyrinth  seal  (p /p  ) 

s  -  labyrinth  seal  intertooth  chamber  width  (in) 

t  -  labyrinth  seal  tooth  thickness  (in) 

t  -  temperature  (°F) 

T-.  -  temperature  upstream  of  orifice  or  nozzle  (°R) 

TQ  -  total  temperature  at  labyrinth  seal  inlet  (°R) 

W  -  mass  flow  rate  (lbm/hr) 

w  -  mass  flow  rate  (lbm/sec) 

w  VTo  |/r7r        ,.         ,  ,, 
w*  =  ~1     ~~- — »     :  non-dimensional  flow  rate 
A  po   f 

X  =  Re  x  10~6 

Y-.  -  expansion  factor,  orifice  or  nozzle 
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a  -  coefficient  of  thermal  expansion,  orifice  or  nozzle  flow  (  F   ) 

0/   -  coefficient  of  discharge,  labyrinth  seal  leakage  (dimensionless) 

6  =  D2/D]_ 

Y  -  ratio  of  specific  heats  =  1.4 

Y*  -  carryover  factor,  labyrinth  seal  leakage 

6  -  labyrinth  seal  tooth  clearance 

A  -  tolerance  designator,  signifies  percent 

co  -  labyrinth  seal  pressure  ratio  function  (dimensionless) 

PI  -  mass  density  (lbm/  ft  ) 

£  -  see  equation  (15) 

Subscripts 

I  -    linear 

n  -  nozzle 

o  -  orifice 

x  -  refers  to  X 

y  -  refers  to  Y 


3.    INSTALLATION: 

The  installation  of  the  Transonic  Turbine  Test  Rig  is  shown  schematically 
in  Figure  1.   Features  of  the  rig  pertinent  to  this  report  are: 

(a)  the  flow  measurement  nozzle  (Figure  2); 

(b)  the  removable  standard  orifice  installation  (see  note  1,  Figure  1 
and  Figure  3) ; 

(c)  the  plenum  labyrinth  leak  rate  measurement  standard  orifice  (2  in. 
pipe  orifice,  Figure  1  and  Figure  4); 

(d)  the  turbine  plenum  labyrinth  seals  (Figure  5) ; 

(e)  the  shaft  labyrinth  seal  (Figure  5) ;  and 

(f)  the  exhauster  installation,  including  associated  piping  (Figure  1). 
This  turbine  can  be  operated  exhausting  either  to  atmospheric  pressure 

or  to  a  partial  vacuum  in  the  hood,  thence  through  the  exhauster.  Figure  1 
shows  all  possible  piping  arrangements,  including  the  removable  calibration 
standard  orifice  installation  and  the  exhauster  with  its  associated  piping. 

The  flow  rate  through  the  turbine  blading  during  operation  is  the  flow 
measured  by  the  flow  nozzle  less  the  plenum  labyrinth  seal  leakage.   Since 
there  is  insufficient  piping  upstream  of  the  nozzle  to  insure  undisturbed  flow, 
the  nozzle  was  calibrated  using  a  standard  orifice  installation,  so  that 
accurate  flow  rates  during  turbine  operation  can  be  obtained  from  the  nozzle 
flow.   Plenum  labyrinth  seal  leak  rates  can  not  be  measured  during  turbine 
operation.   Hence,  such  leakage  is  determined  from  separate  tests  in  which 
all  flow  to  the  plenum  passes  through  the  two  inch  line  standard  orifice 
and  the  only  exit  for  the  flow  is  through  the  labyrinth  seals,  since  the 
turbine  discharge  is  blocked  by  means  of  a  special  cover  plate. 

In  analyzing  exhauster  performance,  the  flow  rates  to  the  exhauster 
from  the  turbine  hood  and  through  the  exhauster  nozzle  are  required.   The 
former  consists  of  the  flow  rates  through  the  turbine,  through  the  labyrinth 


seals  and  through  the  shaft  seals.   There  are  no  provisions  for  actual 
measurement  of  flow  rates  through  the  shaft  seal  or  through  the  exhauster 
nozzle.   However,  with  the  entire  flow  of  the  Allis  Chalmers  compressor 
routed  to  the  turbine  rig,  the  flow  through  the  exhauster  nozzle  is  the 
flow  rate  delivered  by  the  compressor,  determined  as  specified  in  Reference  1, 
minus  the  flow  rate  measured  by  the  flow  nozzle.   In  addition,  since  flow 
through  the  exhauster  nozzle  will  normally  be  choked,  the  flow  rate  can  be 
estimated. 

Locations  of  pressure  taps  and  temperature  probes  are  indicated  in 
Figure  2  for  the  flow  measurement  nozzle,  in  Figure  3  for  the  removable 
standard  orifice  installation,  and  in  Figure  4  for  the  labyrinth  leak  rate 
measurement  orifice.   Both  standard  orifices  are  fitted  with  vena  contracta 
and  flange  taps.   Additionally,  for  both  standard  orifices,  the  lengths  of 
straight  pipe  upstream  and  downstream  of  the  orifices  are  indicated.   Table  I 
lists  applicable  standards  for  obstruction  free,  straight  lengths  of  pipe 
upstream  and  downstream  of  the  orifices  and  for  pressure  tap  locations, 
together  with  the  actual  dimensions  of  the  two  orifices  used  with  the  Transonic 
Turbine  Test  Rig.   Both  installations  are  well  within  prescribed  values  for 
standard  orifice  installations. 

Pressures  were  measured  in  inches  of  mercury  with  atmospheric  pressure 
as  reference;  differential  pressures  were  measured  in  inches  of  water. 
Manometers  used  had  0.1  in  graduations,  with  an  accuracy  of  +  0.03  inches. 
The  accuracy  of  temperature  measurements  was  taken  as  +  1.0  deg  Rankine. 

4.    PROCEDURE: 


Flow  nozzle  calibration  runs  were  made  at  essentially  constant  supply 
pressures  of  either  30  psia  or  45  psia.  Flow  rate  was  controlled  by  using 
the  valve  at  the  exit  of  the  removable  standard  orifice  installation  to  set 


the  differential  pressure  across  the  vena  contracta  taps;  supply  pressure 
was  maintained  by  adjusting  the  Allis  Chalmers  compressor  installation 
main  discharge  (butterfly)  valve.   Differential  pressures  were  increased 
in  increments  of  ten  inches  of  water  for  increasing  flow  rates  to  the 
maximum  differential  pressure,  then  decreased  an  increment  of  five  inches 
of  water,  followed  by  incremental  reductions  of  ten  inches  of  water  to 
minimum  flow.   Approximate  ranges  of  flow  rates  of  the  tests  were  1.44  to 
3.67  lbm/sec  at  30  psia  and  1.8  to  4.6  lbm/sec  at  45  psia. 

Plenum  labyrinth  leak  tests  were  made  by  varying  the  supply  pressure 
from  about  24  to  about  44  psia.   Pressures  were  set  by  observing  the 
turbine  plenum  pressure  reading.   The  pressure  was  varied  by  about  five 
inches  of  mercury  increments  on  pressure  increase,  decreasing  from  maximum 
by  about  2.5  inches  of  mercury  and  then  incremental  decreases  of  about  five 
inches  of  mercury.   Tests  were  made  without  the  exhauster  and  associated 
piping  installed.   Hence  the  maximum  ratio  of  plenum  supply  pressure  and 
labyrinth  discharge  pressure  was  limited  to  approximately  3. 

5.    FORMULAS : 

Square-Edged  Orifice  and  Flow  Nozzle 

The  basic  flow  equation  for  a  square-edged  orifice  or  flow  nozzle, 

in  lbm/hr,  is,  from  page  6,  Reference  2: 

2 


where: 


W  =   359.1   Djf   aKYx    \[p^w  (1) 


D„  =  orifice  diameter  or  nozzle  throat  diameter  (inches) 

2 
a  =  coefficient  of  thermal  expansion  based  on  D   (dimensionless) 
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K  =-'jaischarge  coefficient  (dimensionless) 
Yi  =  expansion  coefficient  (dimensionless) 

P-.  =  fluid  density  (lbm/ft  )  at  upstream  tap 
h  =  differential  pressure  across  taps  (in  H90) 
The  equation  of  state  for  an  ideal  gas  is: 


RT-, 


(2) 


where  pi  is  in  psia,  R  =  53.35  ft  lb/lbm  R,  and  T   is  in  R.   For  P   in  in. 
Hg,  equation  (2)  becomes: 

.  tjjoAna)  (W    .  1328  ;i  (2a) 

(53.35)  Tl  l 

Inserting  (2a)  in  (1),  and  converting  the  latter  to  lbm/sec: 

359.1        2  /  Plhw 

w  ■   3600"     D2   aKYl    J1'328 (3) 

v      T^ 

or: 


Plhw 


w  =  0.115  D2  ciKYjJ  — —  (3a) 

The  coefficient,  a,  is  determined  from  the  linear  coefficient  of  thermal 
expansion  (ot^)  .   For  Do  =  orifice  diameter  at  elevated  temperature: 

'2 
(_2|!=  fl  +  CHtAtYp  =  1  +  2  a.  (At)  +  ap   (At)2     (4) 

1  2        2 
Neglecting  the  higher  order  final  term  and  letting  (D2)   =  o^Dn)  : 

a(D2)2  =  [l  +  2  a£  (At)]  D2  (5) 

and 

a  =  1  +  2  a^  (At)  (6) 

For  both  standard  orifices, made  of  type  304  stainless  steel, 
a£=  (9.5)  10"6  °F-1,  which  gives: 

aQ  =  1  +  0.0019  <jf*j$)  (7) 


(The  valueO. 00193  viceO.0019  was  used  in  calculations,  the  former  value 
being  determined  from  Figure  98,  p.  257,  Ref.  2).   For  the  flow  nozzle 
made  of  2024-T4  aluminum,  o^  =  (12.6)  10~6  °F~1,  and 

an  =  1  +  0.00252  £=&-)  (8) 

The  discharge  coefficient,  K,  is  dependent  on  the  ratio  of  orifice 

or  nozzle  diameter  to  upstream  pipe  diameter  (i.e.,  8  =  D-/D-,)  Reynolds 

Number  based  on  D^,  and  the  type  of  measuring  taps  used.   To  simplify 

computations  for  standard  orifices,  the  following  relation  may  be  used 

(see  p.  65,  Ref.  2) 

K       =       K*  (9) 

[1+A/Re]        [1+A/Re*] 

where,  from  page  212,  Ref.  2: 

A  =  D2  (830-50006+900032-420063+B)  (10) 

with: 

B  =  530/  /—        for  flange  taps  (11) 


B  =  530/,/lK-lOO    for  vena  contracta  taps       (11a) 


and,  from  p.  64,  Ref.  2: 

_6.316W 
e  D2Z*  (12) 

where  Z*  is  absolute  viscosity  in  centipoises  at  flow  temperature  and 
pressure.   For  Z   independent  of  pressure  and  assuming  a  linear  variation 
of  Z   in  the  temperature  range  from  50°F  to  300  F,  there  is  from  p.  335, 
Ref.  2: 

Z=100Z'  =  1.9  +0.24  (ruu  -i  )  (13) 

— f> 
Equation  (13)  in  equation  (12),  with  W=3600w,  and  defining  X=Re  xlO   ,  yields 
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/  -6, 
(6.316)  (3600)  w  (100)   (10  )=   2.275w       (14) 

D2Z  D2Z 


Letting,  Q   -  1  +  A/ 


Re; 

r  -  i  4-  A  10"6 

C  -  X  +  — (15) 


Equation  (9)  becomes: 

K  =  K* 


For  flange  taps,  values  of  K*  =  Koo,  corresponding  to  Re  =°°,  are 
tabulated  in  Table  19  of  Ref.  2.   Equation  (16)  is  then: 

K  =^  Koo  (16a) 

For  vena  contracta  taps,  this  shortcut  method  is  not  authorized  and 
values  of  K  are  not  tabulated.   With  K  tabulated  as  a  function  of  Re  in 

oo 

Table  20  of  Ref.  2,  an  iterative  procedure  using  these  variables  is  normally 
necessary,  with  flow  rate,  w,  as  a  third  variable.   However,  by  judicious 
choice  of  Re   and,  hence,  K  ,  equation  16  becomes  for  vena  contracta  taps: 

K  =  _£_  K*  (16b) 

For  the  removable  standard  orifice  installation,  Re*  was  chosen  as 
(lO  ),  giving  K*=0.698.   For  a  Reynolds  Number  of  2(105),  the  value  of  K 
determined  from  equation  (16b)  differs  by  0.114  percent  from  the  tabulated 
value,  while  the  range  of  Reynolds  Numbers  for  tests  was  approximately: 
3.9  (105)  to  1.3  (106). 

Similarly,  for  the  two  inch  pipe,  square-edged  orifice,  Re*  was  selected 
as  (lO  J,  with  K*  =  0.6103.   At  a  Reynolds  Number  of  5(10^)  the  computed  value 
of  K  differs  from  the  tabulated  value  by  about  0.1  percent  and  at  a  Reynolds 

Number  of  2(10  )  the  difference  is  about  0.13  percent.   The  range  of  Reynolds 

4  5 

Numbers  in  the  tests  conducted  to  date  was  approximately  5(10  )  to  1.2(10  ). 
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Tests  yet  to  be  conducted  with  the  exhauster  installed  will  be  at  somewhat 
higher  Reynolds  Numbers  than  1.2(10  ).   However,  since  it  is  anticipated 
that  the  labyrinth  seal  flow  will  become  choked  at  pressure  ratios  (plenum 
total  to  seal  discharge)  only  slightly  higher  than  3:1,  it  is  considered 
that  Re=2(10  )  is  a  realistic  upper  limit. 

The  expansion  factor,  Yj,  which  compensates  for  compressibility  in  the 
flow  through  the  orifice,  is  defined  by  equation  (60),  p.  71,  Ref.  2,  as: 

Y  =1-(0.41  +  0.3564)  "^ —  (17) 
1                  YPX 

where  Ap  and  p1  must  be  in  consistent  units,  and  y   is  the  specific  heat 
ratio.   With  y   =  1.4,  Ap  =  hw  (in  H-0)  and  p-^   =  P,  (in  Hg) ,  equation  (17) 

becomes: 

Y  =1  Q_  (0.41  ±   0.353*)  h^_  (17a) 
Yl  1'°   (1.4)  (13.59)    v  u/aj 

1 

Calculations  required  to  put  equation  (3a)  in  the  form: 

w  ■  c^p^L  (18) 

for  each  standard  orifice  were  accomplished  by  using  the  CDC  1604  computer. 
Additionally  the  several  constants  in  the  expressions  for  a,  Y.,  and  X  were 
computed.   Table  II  lists  the  values  to  be  used  for  each  standard  orifice. 

It  is  to  be  noted  that  equation  (16a)  was  used  for  flange  taps,  and 
equations  (16b),  (14)  and  (15)  were  used  for  vena  contracta  taps,  along  with 
appropriate  values  of  K^,  K  and  Re  .   In  solving  equation  (18),  £=1  is 
assumed  and  w  is  computed.   With  this  value  of  w,  r  is  obtained  from  equations 
(14)  and  (15),  and  equation  (18)  is  again  solved  using  the  corrected  value  of 

c. 

The  object  of  the  nozzle  calibration  was  the  determination  of  the  discharge 

coefficient  K.   Hence,  equation  (3a)  is  rewritten  as: 

„       w  \Tf7 

Y  V    l     2 ,  (19) 

0.115  D2   anYx  fF~h7 
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where  w  is  the  flow  rate  as  measured  by  the  standard  orifice,  o^  is  given 
by  equation  (8)  and  Y  by  equation  (17a). 

Shaft  Labyrinth  Seal  Leak  Rate 

The  basic  method  for  estimating  the  shaft  labyrinth  seal  leak  rate  is  that 
developed  by  Egli  in  Ref.  3.   The  equation  of  continuity  for  an  ideal  labyrinth 
seal  (i.e.,  all  kinetic  energy  destroyed  after  each  throttling)  is,  from 
Ref.  3: 


A  Po 


w  =  a  /To   V^T 


co 


(20) 


R 


where: 


l-(p/p  ) 

CD  =1       P° 


n  +(2/Y)ln(p0/  ) 
7P 


1/2 


11/2 


1-r' 
T 
n+-  ln(l/r) 


(21) 


a  discharge  coefficient  for  a  single  throttling  (dimensionless) 

2 
A=  cross-sectional  area  of  seal  passage  (in  ) 

p  =  total  pressure  at  seal  entrance  (psia) 

p  =  static  pressure  at  seal  discharge  (psia) 

R=  gas  constant  for  air  (ft  lb/lbm°R) 

To=  total  temperature  at  seal  inlet  (°R) 

(Total  temperature  is  constant  across  seal) 

n=  Number  of  throttlings 

r=  overall  seal  pressure  ratio  (p  /   ) 

po 

The  dimensionless  seal  pressure  ratio  function,  m,  was  developed  by 
series  expansion  with  a  lower  limit  of  pressure  ratio  across  any  one  throttling 
of  0.8.   Since  a  preliminary  analysis  of  exhauster  performance  showed  an  overall 
seal  pressure  ratio  of  about  0.55  (hood  case  pressure  to  ambient  pressure),  this 
same  assumption  was  made  for  the  ten  throttlings  through  the  shaft  labyrinth  seal. 

To  account  for  a  non-ideal  labyrinth,  a  carryover  factor,  y    ,  is  introduced 
to  adjust  for  carryover  of  kinetic  energy  from  one  throttling  to  the  next. 
Equation  (20)  then  becomes: 
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w=ay  -^p — Jg/    co  (22) 


where 

a=f(S/t,t)  =  constant 

y*  =  f(5/s»  n)  =  constant 

Here,  6  =  tooth  clearance,  t=tooth  width  and  s=chamber  width  (see  Fig.  5). 

Unfortunately,  the  curves  of  a  versus  6/t  given  in  Ref.  3  make 
determination  of  a  for  the  geometry  of  the  shaft  seal  somewhat  uncertain. 
Hence,  the  value  of  a=0.76  was  taken  from  experimental  data  given  on  page  72, 
Ref.  4,  where  a  is  plotted  as  a  function  of  6/t  and6/b,  where  b  is  the  chamber 
heighth. 

Plots  of  y*   vs  6/s  given  in  Ref.  3  are  based  on  experimental  data  for 
6>  0.01.   As  cited  in  that  reference,  Friedrich's  tests  of  straight  through 
labyrinths  with  6=  0.006  to  0.010  inches,  and  with  two  and  three  throttlings, 
indicated  somewhat   smaller  values  of  y*   than  those  given  in  Ref.  3  for  the 
geometry  of  the  shaft  seal  of  the  Transonic  Turbine  Test  Rig.   Considering, 
also,  ten  rather  than  two  or  three  throttlings  and  data  for  multiple  throttlings 
given  in  Ref.  4,  Friedrich's  value  of  y*   =  1.15  was  selected  in  preference  to 
Egli's  value  of  y*=   1.25. 

With  these  values  of  a  and  y* ,    the  final  formula  for  estimating  shaft 
seal  leakage  is,  with:   A=   n  1.2  x  .005: 


w=  0.0128  ^r-   <p  (23) 


With  w*  defined  by 


w*  =  w  U°       VR/g  =   0.874  cp  (24) 

A  P  o 


=  ^E-^R/g  =   0.874  cp 
A  P  o 

a  single  plot  of  w*  versus  overall  seal  pressure  ratio  (r)  can  be  used  since: 

cp  =  f(r) 

Plenum  Labyrinth  Seal  Leak  Rate 

Since  plenum  labyrinth  leak  rates  were  determined  experimentally,  use  of 
equation  (22)  was  limited  to  attempts  to  correlate  experimental  data  with  the 
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method  and  to  extrapolate  higher  overall  seal  pressure  ratio  data  (0.34  to 
0.6)  to  lower  pressure  ratios  (0.16  to  0.34).   For  the  two  identical  labyrinths 
on  either  side  of  the  plenum,  equation  (22)  is  written: 

w*  =  "^   |/r7=  2Kcp  (25) 

A  -o     '   g 

where  K=ay*  of  equation  (22).   Then,  if  a  and  y*  (or  ay*)  are  in  fact  constant, 

K=constant,  and  experimental  data  at  higher  values  of  r  can  be  extrapolated 

to  lower  values  of  r.   Since  it  is  anticipated  that  flow  through  the  seal  will 

be  choked  at  (/max  =  0.2798  (r=  0.225),  the  range  of  extrapolation  would  be 

relatively  small. 

6.    RESULTS: 


Flow  Nozzle  Calibration 

Since  the  tolerance  for  vena  contracta  taps  is  less  than  that  for  flange 
tops  (see  Section  8),  plotted  and  tabulated  data  are  vena  contracta  values 
unless  otherwise  indicated. 

Figure  6  shows  plots  of  nozzle  discharge  coefficient,  Kn,  versus  Reynolds 
Number  determined  from  tests  of  8  September  and  13  September  1965.   Figure  6 
indicates  that  Kn  can  be  taken  as  a  constant  above  Reynold's  Numbers  of  about 
6(10  )  -  (w=21bm/sec) .   By  the  method  of  least  squares,  values  of  Kn=constant 
were  determined  for  the  following  data  groups:   8  September  tests  at  30  psia; 
13  September  tests  at  30  psia;  these  two  tests  combined;  13  September  tests  at 
45  psia;  and  all  the  above-listed  tests  of  8  and  13  September  combined.   Re- 
sulting minimum  error  values  of  Kn=constant  are  tabulated  in  Table  III  with 
the  least  squares  error.   In  addition,  certain  values  of  Kn=constant  are  shown 
in  Figure  6.   Values  in  Table  III  show  reasonably  close  agreement.   However, 
the  13  September  test,  vena  contracta  values  are  considered  preferable,  since 
vena  contracta  flow  rates  are  generally  more  accurate  than  flange  tap  flow 
rates,  and  since  some  difficulty  was  encountered  with  water  collecting  in 
orifice  pressure  tap  lines  during  the  8  September  tests.   (Even  though  these 
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lines  were  continually  drained  during  the  8  September  tests,  the  data  are 
considered  inferior  to  those  from  the  13  September  tests,  during  which 
the  water  collection  problem  did  not  exist.) 

Shaft  Labyrinth  Seal  Leak  Rate 

Figure  7  shows  w*  and  (/^plotted  versus  overall  seal  pressure  ratio  (r)- 
see  equation  (21)  -  for  use  in  estimating  shaft  seal  leakage. 

Plenum  Labyrinth  Seal  Leak  Rate 

Experimentally  determined  values  of  K  and  w*  from  equation  (25)  are  plotted 
versus  overall  seal  pressure  ratio  in  Figures  8  and  9  respectively.   It  is 
apparent  from  Figure  8  that  extrapolation  of  K  to  lower  pressure  ratios  would 
be  of  questionable  accuracy.   To  provide  a  mathematical  expression  for 
w*=f(r),  the  data  shown  in  Figure  9  was  fit  with  a  parabola  by  the  method  of 
least  squares.   In  this  curve  fitting,  a  parabola  of  the  form: 

r=  a+bw*+cw*  (26) 

was  assumed.   In  addition  to  experimental  data,  the  known  condition  that  at 
r=1.0,  w*=0.0,  was  introduced.   Symmetry  about  the  r  axis  was  assumed  by 
including  in  the  least  squares  formulation  both  plus  and  minus  values  of  w* 
for  each  data  point  used  (i.e.,  b=0.0). 

During  the  26  October  tests,  especial  care  was  taken  to  assure  flow 
stabilization  before  data  was  taken.   Although  the  difference  between  increasing 
pressure  and  decreasing  pressure  data  was  not  eliminated,  it  was  materially 
reduced  from  this  difference  for  the  20  October  tests.   Consequently,  the 
decreasing  pressure  data  (points  to  the  far  right  in  Figures  8  and  9)  were 
neglected  in  establishing  both  the  solid  curves  of  Figures  8  and  9  and  the 
least  squares  parabola. 

The  resulting  least  squares  second  order  polynomial,  for  which  the  least 
squares  is  0.02079,  is: 


or,    solved   for  w 


* 


r=1.0-4.0791w*2  (27) 


Jl.O-r 
w*  =    1/4.0791  (27a) 


This  equation  is  plotted  in  Figures  9  and  10  along  with  experimental 
data.  Also  shown  in  Figure  10  are  segments  of  three  curves,  w*=2K^, 
for  values  of  K  appropriate  to  three  different  ranges  of  seal  pressure  ratio. 

7.    DISCUSSION  AND  RECOMMENDATIONS: 

Based  on  the  nozzle  flow  rate  experiments  conducted,  it  is  recommended 
that  a  valve  of  0.9992  be  used  for  the  nozzle  discharge  coefficient,  Kn,  in 
equation  (3a).   However,  when  supply  pressures  to  the  nozzle  are  30  psia  or 
below  a  valve  of  1.0012  should  be  used. 

Equation  (27a)  is  recommended  for  determination  of  plenum  labyrinth  leak 
rates  in  the  pressure  ratio  range:   0.03  to  1.0.   Since  plenum  labyrinth  leak 
rates  (0.04  to  0.10  lbm/sec)  are  very  much  smaller  than  turbine  flow  rates 
(2.0  to  3.9  lbm/sec),  the  slight  errors  in  the  leakage  flow  rates  caused  by 
using  equation  (27a)  should  have  a  negligible  effect  on  the  accuracy  of  turbine 
flow  rates.   Furthermore,  this  equation  should  give  satisfactory  flow  rates  at 
higher  pressure  ratios  (see  Figure  10). 

For  exhauster  operation  with  a  supply  pressure  of  45  psia  to  the  exhauster 
nozzle  and  the  turbine,  the  estimated  shaft  seal  leak  rate  given  by  equation 
(24)  is  about  0.003  lbm/sec,  while  flow  rates  through  the  exhauster  nozzle  and 
the  turbine  are  estimated  as  9.38  and  3.87  lbm/sec  respectively.   Using  Egli's 
values  of  nozzle  discharge  coefficient  and  carryover  factor  (a==  0.9  and  y*=1«25), 
shaft  seal  leak  rate  estimates  would  be  approximately  50%  greater  than  estimates 
obtained  using  equation  (24).   However,  plenum  labyrinth  leak  rates  predicted 
using  Egli's  values  of  a  and  y*  are  about  50  to  70  percent  greater  than  the 
leak  rates  determined  experimentally.   Hence,  the  use  of  equation  (24)  is 
recommended  for  predicting  shaft  seal  leakage  flow  rates. 
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The  recommended  flow  rate  equations  are  summarized  in  Table  IV  for 
ready  reference. 


8.         TOLERANCE   IN   FLOW  MEAgTJRI 


*— «— 


4ymmA#ta*» 


The  overall  error  tolerance  in  percent  in  flow  measurements  using 
sharp-edged  orifices  is  given  by  the  square  root  of  the  sum  of  the  squares 
of  "exponents  times  tolerances,  in  percent"  ascribable  to  the  individual 
items  in  the  flow  equation.   Hence,  from  page  24,  Ref.  2: 

A(w)  =  |[2A(D2)]  2  +  [A(c.)]  2  +  [a«0]  2  +  [a(Yx)]  2 

+  f|A(hn)]  2+[^A(P1)]  2  +  [^MT1)]2)|  (28) 

Tolerances  of  individual  items  for  the  Transonic  Turbine  Test  Rig  installation, 
based  on  Reference  2,  are  given  in  Table  V. 

The  error  in  values  of  discharge  coefficients  of  the  flow  nozzle  and  the 
sharp-edged  orifice  must  be  approximately  the  same.   It  follows  also  that  the 
tolerances   in  flow  nozzle  flow  rate  measurements  must  be  the  same  as  those  for 
the  sharp-edged  orifice  used  in  calibration.   Hence,  vena  contracta  tap  based 
values  are  preferable  due  to  a  lower  tolerance  than  that  for  flange  taps. 

Tolerances  in  measured  labyrinth  leakage  rates  also  should  be  in  accordance 
with  the  preceding  paragraph.   However,  it  is  apparent  from  Figures  8  and  9  that 
there  is  considerable  scatter  in  the  data.   Hence,  a  tolerance  arbitrarily  set 
at  +  10%  would  appear  reasonable. 

Tolerances  in  estimated  shaft  seal  leak  rates  are  difficult  to  specify 
exactly.   Here,  a  reasonable  tolerance  could  be  set  at  +  20%,  based  on 
qualifications  set  forth  in  the  literature  cited. 
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TABLE  II 


Flow  Equation  Constants  -  Square-Edged  Orifices 
Transonic  Turbine  Test  Rig 


Orifice 
Installation 

6 

Flange 

Taps 

Vena  Cont  Taos 

Cx 

._ — i 

Cy 

C 

A 

c 

A 

EQ(18) 

EQ(10) 

EQ(18) 

EQ(10) 

EQ(14) 

EQ(17a) 

Removable 

0.6995 

1.4257 

0.00218 
y  itfi 

1.4405 

0.001756 

.5359 
| 

0.026 

2"  pipe 

0.3991 

0.0476 

0.000305 

0.04761 

0.0002249 



:  2.7561 

I  0.022 

EQUATION 
(18) 

(10) 
(14) 

(17a) 

(13) 

(8) 
(15) 


£«§*! 


Pl  hw 


A  =  D  (830-50006  +  90006  -  42006  +  B) 
-4 


=  (6.316)(3600U0  w   =  c  w 


D2Z 


T  -  1.0- 


(.41  +  .356  ) 


1.4  x  13.59 


x  Z 


h  hw 

w    =  x_c  — 

-^~  y   Pi 


■  1-9  +  0-24(rdo-1} 

\  100/ 


a  =  1.0  +  0.00252 


C   =  1.0  + 


a(io-^) 


1 


TABLE  III 

Flow  Nozzle  Discharge  Coefficients 
Transonic  Turbine  Test  Rig 


Test 

Nominal 
Pressure 

Data 
Points 

Flange 

Taps 

Vena  Cont  Taps 

Date 

Error  # 

Kn 

Error  # 

Kn 

8  Sep 

30  psia 

9 

.00828 

.9994 

.00599 

1.0038 

13  Sep 

30  psia 

12 

.01055 

1.0026 

.00684 

1.0012* 

8/13  Sep 

30  psia 

21 

.01531 

1.0012 

.01071 

1.0024 

8/13  Sep 

30  psia 
45  psia 

21 
12 

.01614 

1.0011 

.01633 

1.0012 

13  Sep 

1 
45  psia 

12 

.00497 

1.0008 

-  —  i 
.00864 

.9992* 

//   Error  figures  show  least  squares  average  difference  of  data  points  from 
Kn  for  minimum  error. 

*  Recommended  values 
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TABLE  IV 

Summary  of  Formulas  for  Determining  Flow  Rates 
Transonic  Turbine  Flow  Rates 


Desired 
Flow  Rate 


Basic 
Formula 


Auxiliary 
Formulas 


Nozzle 
(lbm/sec) 


w=2.073anY1 


Pl  hw 


a=l. 0+0. 00252 


t1-68 


100 


Y]  =  1.0-0.023   —" 


Plenum  Labyrinth 
Seal  Leak  Rate 
(lbm/sec) 


w=0.116    po 


fTc 


I  1.0-r 
w*  =V   4.07< 


0791 


r  = 


Po 


Shaft  Labyrinth 
Seal  Leak  Rate 
(lbm/sec) 


w  -   0.01466 


w*po 


w*   =   0.874  <P 
.2 


1/2 


V- 


1.0-r' 


}°'+—h(l/r). 


po 


r 


cp_ 


25 


TABLE  V 

Percent  Tolerances  of  Flow  Equation  Variables 
and  Parameters  -  Square-Edged  Orifices  * 

Transonic  Turbine  Test  Rig 


a)   Independent  of  Orifice: 
A  (a)  =  +  .05 
A  (Yp  =  0  to  +  .05 

=  +  0.5  to  +  1.0 


A  ,,  v    +  .03 

A  (hw)  =  - (ioo) 


hw 


A  (P]_)  =  -   *03    (100) 
~1 

A  (j ^  m  ±  i-o   (100) 


hw 


for  0.0  <  PY  <  -015 
for  0.015  <  Yi 


b)   Removable  Standard  Orifice  Installation: 
0005 


A  (D2}  =±  4.2425 


(100)  =  +  .0118 


A  (C)  =  +  1.0 

=  +  0.4  to  +  0.6 
c)   2"  Pipe  Orifice 

A  (D2)  =  ±7§25^  (100>  "  -0606 

A  (C)  =  +  1.0 

=  +  0.4  to  +  0.6 


flange  taps 

vena  contracta  taps 


flange  taps 

vena  contracta  taps 


*  Basis:   Ref.  2 
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U.  S.  Naval  Postgraduate  School 


Piping  Installation 
Transonic  Turbine  Test  Rig 


FIG.    1 


TECH.    NOTE" 
N0.     66T-1 


Orifice 


Flow 
Siraiyhtener 


tank  Z 


L-Uj    / —  Exhauster  (fche,  4) 


Shaft   Seal 


NOTES 

1.  Installed   for    Nozzle 
Calibration    Only. 

2.  Blocked  During  Nozzle 
Calibration.. 

3.  Blocked   During    Laby- 
rinth  Leak  Tests. 

4.  Exhauster   and  Associ- 
ated   Piping    Required 
Only    fcr    Pressure    Ra- 
tios  Greater  Than, 3'- 1. 


Inlet 
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Flow  Measurement  Nozzle 
Transonic  Test  Rig 


FIG.  2 
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Removable  Standard  Orifice  Installation 
Transonic  Turbine-  Test  Rig 


FIG.  3 


TECH  NOTE 
N0.66T-1 
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2"  PIPE  STANDARD  ORIFICE  INSTALLATION 
Transonic  Turbine  Test  Rig 


FIG.    4 
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